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Abstract 

Superlattice provides a new approach to enrich the class of materials with novel properties. Here, we report the 
structural and electronic properties of superlattices made with alternate stacking of two-dimensional hexagonal 
germanene (or silicene) and a MoS 2 monolayer using the first principles approach. The results are compared with 
those of graphene/MoS 2 superlattice. The distortions of the geometry of germanene, silicene, and MoS 2 layers due 
to the formation of the superlattices are all relatively small, resulting from the relatively weak interactions between 
the stacking layers. Our results show that both the germanene/MoS 2 and silicene/MoS 2 superlattices are manifestly 
metallic, with the linear bands around the Dirac points of the pristine germanene and silicene seem to be preserved. 
However, small band gaps are opened up at the Dirac points for both the superlattices due to the symmetry breaking in 
the germanene and silicene layers caused by the introduction of the MoS 2 sheets. Moreover, charge transfer happened 
mainly within the germanene (or silicene) and the MoS 2 layers (intra-layer transfer), as well as some part of the 
intermediate regions between the germanene (or silicene) and the MoS 2 layers (inter-layer transfer), suggesting 
more than just the van der Waals interactions between the stacking sheets in the superlattices. 
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Background 

In the past decade, the hybrid systems consisting of gra- 
phene and various two-dimensional (2D) materials have 
been studied extensively both experimentally and theor- 
etically [1-6]. It has long been known that the thermal, 
optical, and electrical transport properties of graphene- 
based hybrids usually exhibit significant deviations from 
their bulk counterparts, resulting from the combination 
of controlled variations in the composition and thickness 
of the layers [6,7]. Moreover, the use of 2D materials 
could be advantageous for a wide range of applications 
in nanotechnology [8-13] and memory technology 
[14-16]. Among those hybrid systems, the superlattices 
are considered as one of the most promising nanoscale 
engineered material systems for their possible applications 
in fields such as high figure of merit thermoelectrics, 
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microelectronics, and optoelectronics [17-19]. While the 
research interest in graphene-based superlattices is growing 
rapidly, people have started to question whether the gra- 
phene could be replaced by its close relatives, such as 2D 
hexagonal crystals of Si and Ge, so called silicene and ger- 
manene, respectively. Silicene and germanene are also zero- 
gap semiconductors with massless fermion charge carriers 
since their n and if bands are also linear at the Fermi level 
[20]. Systems involving silicene and germanene may also be 
very important for their possible use in future nanoelectro- 
nic devices, since the integration of germanene and silicene 
into current Si-based nanoelectronics would be more likely 
favored over graphene, which is vulnerable to perturbations 
from its supporting substrate, owing to its one-atom 
thickness. 

Germanene (or silicene), the counterpart of graphene, 
is predicted to have a geometry with low-buckled honey- 
comb structure for its most stable structures unlike the 
planar one of graphene [20-22]. The similarity among 
germanene, silicene, and graphene arises from the fact 
that Ge, Si, and C belong to the same group in the 
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periodic table of elements, that is, they have similar elec- 
tronic configurations. However, Ge and Si have larger 
ionic radius, which promotes sp 3 hybridization, while 
sp 2 hybridization is energetically more favorable for C 
atoms. As a result, in 2D atomic layers of Si and Ge 
atoms, the bonding is formed by mixed sp 2 and sp 3 
hybridization. Therefore, the stable germanene and sili- 
cene are slightly buckled, with one of the two sublattices 
of the honeycomb lattice being displaced vertically with 
respect to the other. In fact, interesting studies have 
already been performed in the superlattices with the in- 
volvement of germanium or/and silicon layers recently. 
For example, the thermal conductivities of Si/SiGe and 
Si/Ge superlattice systems are studied [23-25], showing 
that either in the cross- or in-plane directions, the sys- 
tems exhibit reduced thermal conductivities compared 
to the bulk phases of the layer constituents, which im- 
proved the performance of thermoelectric device. It is 
also found that in the ZnSe/Si and ZnSe/Ge superlat- 
tices [26], the fundamental energy gaps increase with 
the decreasing superlattice period and that the silicon 
or/and germanium layer plays an important role in de- 
termining the fundamental energy gap of the superlat- 
tices due to the spatial quantum confinement effect. 
Hence, the studies of these hybrid materials should be 
important for designing promising nanotechnology 
devices. 

In the present work, the structural and electronic 
properties of superlattices made with alternate stacking 
of germanene and silicene layers with MoS 2 monolayer 
(labeled as Ger/MoS 2 and Sil/MoS 2 , respectively) are 
systematically investigated by using a density functional 
theory calculation with the van der Waals (vdW) correc- 
tion. In addition, we compare the results of Ger/MoS 2 and 
Sil/MoS 2 superlattices with the graphene/MoS 2 superlattice 
[6] to understand the properties concerning the chemical 
trend with the group IV atoms C, Si, and Ge in the super- 
lattices. Our results show that Ger/MoS 2 and Sil/MoS 2 
consist of conducting germanene and silicene layers and 
almost-insulating MoS 2 layers. Moreover, small band gaps 
open up at the K point of the Brillouin zone (BZ), due 
to the symmetry breaking of the germanene and sili- 
cene layers which is caused by the introduction of the 
MoS 2 layers. Localized charge distributions emerged 
between Ge-Ge or Si-Si atoms and their nearest neigh- 
boring S atoms, which is different from the graphene/ 
MoS 2 superlattice, where a small amount of charge 
transfers from the graphene layer to the MoS 2 sheet 
[6]. The contour plots for the charge redistributions 
suggest that the charge transfer between some parts of 
the intermediate regions between the germanene/sili- 
cene and the MoS 2 layers is obvious, suggesting much 
more than just the van der Waals interactions between 
the stacking sheets in the superlattices. 



Methods 

The present calculations are based on the density func- 
tional theory (DFT) and the projector-augmented wave 
(PAW) representations [27] as implemented in the 
Vienna Ab Initio Simulation Package (VASP) [28,29]. 
The exchange-correlation interaction is treated with the 
generalized gradient approximation (GGA) which is pa- 
rameterized by Perdew-Burke-Ernzerhof formula (PBE) 
[30]. The standard DFT, where local or semilocal func- 
tionals lack the necessary ingredients to describe the 
nonlocal effects, has shown to dramatically underesti- 
mate the band gaps of various systems. In order to have 
a better description of the band gap, corrections should 
be added to the current DFT approximations [31,32]. 
On the other hand, as is well known, the popular dens- 
ity functionals are unable to describe correctly the vdW 
interactions resulting from dynamical correlations be- 
tween fluctuating charge distributions [33]. Thus, to im- 
prove the description of the van der Waals interactions 
which might play an important role in the present lay- 
ered superlattices, we included the vdW correction to 
the GGA calculations by using the PBE-D2 method [34]. 
The wave functions are expanded in plane waves up to a 
kinetic energy cutoff of 420 eV. Brillouin zone integrations 
are approximated by using the special /c-point sampling of 
Monkhorst-Pack scheme [35] with a T-centered 5x5x3 
grid. The cell parameters and the atomic coordinates of the 
superlattice models are fully relaxed until the force on each 
atom is less than 0.01 eV/A. 

Results and discussions 

For the free-standing low-buckled germanene and sili- 
cene, the calculated lattice constants are 4.013 and 
3.847 A, respectively, which agree well with the reported 
values of 4.061 and 3.867 A for germanene and silicene, 
respectively [36]. Our optimized lattice constant for a 
MoS 2 monolayer is 3.188 A, which is the same as the 
previous calculated values by PBE calculations [37]. Al- 
though the lattice constants of germanene/silicene and 
MoS 2 monolayer are quite different, all of them do share 
the same primitive cell of hexagonal structure. For es- 
tablishing the calculation models for Ger/MoS 2 and Sil/ 
MoS 2 superlattices and to minimize the lattice mismatch 
between the stacking sheets, we have employed super- 
cells consisting of 4 x 4 unit cells of germanene (and sili- 
cene) and 5x5 unit cells of MoS 2 monolayer in the x-y 
plane. Thus, we have 4^(gemanene) = 16.052 A, 4^(sili- 
cene) = 15.388 A, and 5^(MoS 2 monolayer) = 15.940 A, 
which lead to a lattice mismatch of around 0.70% 
between the germanene and MoS 2 layers and 3.46% be- 
tween the silicene and MoS 2 layers. Compared with the hy- 
brid systems investigated previously [38-42], the present 
lattice mismatch values are very small. In the calculations, 
first, the lattice constant of germanene/silicene (4<2 ger/sil ) 
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was set to match to that (5^moS 2 ) of the MoS 2 monolayer 
in the supercell The supercells are then fully relaxed for 
both the lattice constants and the atomic geometry. The 
mismatch will finally disappear, leading to the commensur- 
ate systems. The superlattices we introduced in this work, 
by hybridizing germanene or silicene with MoS 2 monolayer, 
are shown in Figure 1. The supercells consist of alternate 
stacking of one germanene or silicene sheet and one MoS 2 
monolayer, with 32 Ge or Si atoms, 25 Mo, and 50 S atoms 
per supercell. For a single Ge or Si atom adsorbed on a 
MoS 2 monolayer, there are three possible adsorption sites, 
i.e., the top site directly above a Mo atom, the top site dir- 
ectly above a S atom, and the hollow site above the center 
of a Mo-S hexagon. For the Ger/MoS 2 and Sil/MoS 2 super- 
lattices, we consider two possible representative arrange- 
ments of germanene/silicene on the MoS 2 monolayer: (i) 
one Ge or Si atom in the supercell (4x4 unit cell) was set 
to sit directly on top of one Mo/S atom (the positions of all 
the other Ge or Si atoms will then be determined). In this 
way, there will be one Ge or Si atom in the supercell sitting 
on top of a S/Mo atom, too; see Figure lc. (ii) One Ge or Si 
atom in the supercell was set to sit on the hollow site above 
the center of a hexagon of MoS 2 , as shown in Figure Id. 
From the present calculations, it is found that the binding 
energy differences between the above models of superlat- 
tices are very small (about 1 to 2 meV), which indicates that 
the energy of superlattice is not sensitive to the stacking of 
the atomic layers. Thus, in this paper, we show only the re- 
sults of the configuration with one Ge or Si atom on top of 
the Mo or S atom. In all the stacking types, the 2D 



characteristics of the superlattice structures are kept, e.g., 
hexagonal atomic networks are seen in both Figure lc,d 
which shows the fully optimized geometric structures of 
the supercells. Actually, the changes of the superlattice 
structures are quite small by atomic relaxations. The calcu- 
lated lattice constants of Ger/MoS 2 and Sil/MoS 2 superlat- 
tices are 15.976 and 15.736 A, respectively. In the Ger/ 
MoS 2 superlattice, the germanene layers are compressed by 
0.47% (from 4.013 to 3.994 A) as compared to the corre- 
sponding isolated germanene, while the MoS 2 layers are ex- 
panded by 0.22% (from 3.188 to 3.195 A) as compared to 
the free-standing MoS 2 monolayer. On the other hand, in 
the case of Sil/MoS 2 superlattice, the silicene layers in the 
superlattice are expanded by 2.26% (from 3.847 to 3.934 A), 
while the MoS 2 layers in the supercell are reduced by 1.29% 
(from 3.188 to 3.147 A) (see Table 1). 

The averaged Mo-S bond lengths of the superlattices 
are calculated to be all around 2.400 A (see Table 1). 
The averaged Ge-Ge/Si-Si bond lengths (d Ge _ G J d S isd in 
the relaxed superlattices are all around 2.400/2.300 A, 
which are close to those in the free-standing germanene/ 
silicene sheets (2.422/2.270 A). Although the atomic bond 
lengths in the stacking planes are almost the same for 
Ger/MoS 2 and Sil/MoS 2 superlattices, the interlayer dis- 
tances (d) exhibit relatively larger deviations (but still close 
to each other; see Table 1). A shorter interlayer distance d 
is found in the Ger/MoS 2 system, indicating that the 
Ge-MoS 2 interaction is stronger than the Si-MoS 2 inter- 
action in the Sil/MoS 2 system. The Ge-S and Si-S atomic 
distances in the Ger/MoS 2 and Sil/MoS 2 superlattices are 
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Table 1 Binding energies, geometries, supercell lattice constants, averaged bond lengths, sheet thicknesses, and 
buckling of superlattices 



System 


E b (per Ge/Si) 


E b (per MoS 2 ) 


a = b 


c 


dittos 


^Ge-Ge/dsi-Si 


h s -s 


^Ge 






(eV) 


(eV) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


Ger/MoS 2 


0.277 


0.354 


15.976 


9.778 


2.410 to 2.430 


2.420 to 2.440 


3.129 


0.782 




Sil/MoS 2 


0.195 


0.250 


15.736 


9.926 


2.400 to 2.410 


2.320 to 2.330 


3.176 




0.496 


Germanene 






16.052 






2.422 




0.706 




Silicene 






15.388 






2.270 






0.468 


MoS 2 monolayer 






15.940 




2.413 




3.118 







Theoretical geometries of the isolated germanene, silicene, and MoS 2 monolayer are also listed. E b , binding energies (per Ge/Si atom and per MoS 2 ); a, b, and c, 
supercell lattice constants; c/ Mo . s , d Ge . Ge , and d S \s], averaged Mo-S and Ge-Ge/Si-Si bond lengths; h S s, sheet thicknesses of MoS 2 ; A Ge and Z\ Si , amplitude of buckling 
of the germanene and silicene in the superlattices. 



2.934 and 3.176 A, respectively, where both values are 
shorter than 3.360 A in the graphene/MoS 2 superlattice [6]. 
Such decreases of interlayer distances indicate the enhance- 
ment of interlayer interactions in the Ger/MoS 2 and 
Sil/MoS 2 superlattices as compared to the graphene/MoS 2 
one. This can also explain why the amplitude of buckling 
(A) in the germanene/silicene layers of the superlattices be- 
come larger as compared to the free-standing germanene/ 
silicene, i.e., A going from 0.706 to 0.782 A in the germa- 
nene layers and from 0.468 to 0.496 A in the silicene layers. 
The Ge-S and Si-S atomic distances in the Ger/MoS 2 and 
Sil/MoS 2 superlattices (2.934 and 3.176 A) are much larger 
than 2.240 and 2.130 A, the sum of the covalent atomic ra- 
dius of Ge-S and Si-S atoms (the covalent radius is 1.220/ 
1.110 A for germanium/silicon and 1.020 A for sulfur), 
which suggests that the interlayer bonding in the superlat- 
tices is not a covalent one. 

To discuss the relative stabilities of the superlattices, 
the binding energy between the stacking sheets in 
the superlattice is defined as E\> = - [E SU perceii -(^MoS 2 + 
^Ger/Sii)]/-W > where Supercell is the total energy of the 
supercell, and £moS 2 and E Ger/Si \ are the total energies of 
a free-standing MoS 2 monolayer and an isolated germa- 
nene/silicene sheet, respectively. When N = N(Ge/Si) = 
32, the number of Ge/Si atoms in the supercell, E h is 
then the interlayer binding energy per Ge/Si atom. 
When N = Af(MoS 2 ) = 25, the number of sulfur atoms in 
the supercell, then, E h is the interlayer binding energy 
per MoS 2 . The interlayer binding energies per Ge/Si 
atom and those per MoS 2 are presented in Table 1. 
£moS 2 i s calculated by using a 5 x 5 unit cell of the MoS 2 
monolayer, and E Ger/Si i is calculated by using a 4 x 4 unit 
cell of the germanene/silicene. The binding energies be- 
tween the stacking layers of the superlattices, calculated 
by the PBE-D2 method, are both relatively small, i.e., 
0.277 eV/Ge and 0.195 eV/Si for the Ger/MoS 2 and Sil/ 
MoS 2 superlattices, respectively (see Table 1). The small 
interlayer binding energies suggest weak interactions be- 
tween the germanene/silicene and the MoS 2 layers. The 
binding energy also suggests that the interlayer interaction 



in Ger/MoS 2 superlattice is slightly stronger than that in 
the Sil/MoS 2 one. The interlayer binding energies are 
0.354 eV/MoS 2 and 0.250 eV/MoS 2 for the Ger/MoS 2 and 
Sil/MoS 2 superlattices, respectively, both are larger than 
0.158 eV/MoS 2 in the graphene/MoS 2 superlattice [6]. This 
is an indication that the mixed spi-sp-i hybridization in the 
buckled germanene and silicene leads to stronger bindings 
°f germanene/silicene with their neighboring M0S2 atomic 
layers, when compared with the pure planar spi bonding in 
the graphene/MoS2 superlattice. In addition, the interlayer 
bindings become stronger and stronger in the superlattices 
of graphene/MoS2 to silicene/MoS2 and then to germa- 
nene/MoS2 monolayer. 

Figure 2 shows the band structures of various 2D ma- 
terials, e.g., the bands of flat germanene/silicene com- 
pared with low-buckled germanene/silicene. The band 
structure of flat silicene is similar to that of low-buckled 
one. In both kinds of silicene, the systems are semimetal 
with linear bands around the Dirac point at the K point 
of the Brillouin zone. On the other hand, the band struc- 
ture of flat germanene is quite different from that of 
low-buckled one. The flat germanene is metallic, and the 
Dirac point does not sit at the Fermi level (but above 
the E F ). The band structure of low-buckled germanene, 
however, is similar to that of the low-buckled silicene. 
To help understand the electronic band structures of the 
superlattices and the contribution of each atomic layer 
to the band structures, we present in Figure 3 the band 
structures of Ger/MoS 2 and Sil/MoS 2 superlattices, to- 
gether with those of the independent low-buckled ger- 
manene/silicene and MoS 2 monolayer sheets. The band 
structures of free-standing buckled germanene/silicene 
and MoS 2 sheets (Figure 3a,b,c) are calculated by using 
4x4 and 5x5 supercells, respectively, in order to com- 
pare with the band structures of the superlattices dir- 
ectly. The band structures of the Ger/MoS 2 and Sil/ 
MoS 2 superlattices are presented in Figure 3d,e, where 
the contributions of the germanene/silicene and MoS 2 
monolayers to the band structures of the superlattices 
are shown with blue and green dots (where the size of 



Li et a I. Nanoscale Research Letters 2014, 9:1 10 
http://www.nanoscalereslett.eom/content/9/1/1 10 



Page 5 of 9 




dots are proportional to the contributions), respectively. 
In general, the outlines of the band structures of the two 
superlattices seem to be similar to the rigid sum' of the 
bands of each constituent (i.e., the bands of independent 
germanene/silicene and MoS 2 sheets), indicating that the 
couplings between the stacking sheets are relatively 
weak. However, new important characters in the band 
structures of the superlattices appear. Both the Ger/ 
MoS 2 and Sil/MoS 2 superlattice systems manifest metal- 
lic properties, since there are several bands crossing the 
Fermi level. In fact, in the superlattice systems, the Dirac 
points of the free-standing germanene/silicene (at the K 
point) move upward slightly above the Fermi level; at 
the same time, the Dirac points at the H point (H is 
above K in the ^-direction in the BZ) move downward 
slightly below the Fermi level. Such shifts of Dirac points 
lead to partially occupied bands in the superlattices, also 
implying charge transfer around K point to the H point 
in the BZ. The bands crossing the Fermi level are con- 
tributed mainly by the germanene/silicene layers rather 
than the MoS 2 sheets in both the Ger/MoS 2 and Sil/ 
MoS 2 superlattices, except that small contributions from 
MoS 2 sheet are visible around the H point. Contribu- 
tions from the MoS 2 layers to the electronic states 



around the Fermi level are more significantly visible in 
the system of Ger/MoS 2 than in the Sil/MoS 2 system. 
The feature of energy bands suggests that the electronic 
conduction of the superlattices exists mainly in the x-y 
plane and is almost contributed by the germanene/sili- 
cene sheets rather than the MoS 2 sheets, namely, the 
superlattices are compounds made with alternate stack- 
ing of conductive germanene/silicene layers and nearly 
insulating MoS 2 sheets. This is different from the gra- 
phene/MoS 2 superlattice, in which both graphene and 
MoS 2 layers can be conductive, resulting from the 
charge transfer between the graphene and MoS 2 sheets 
[6]. Moreover, according to the detailed band structures 
inserted in the vicinity of Figure 3d,e, we found that 
small band gaps opened up at the K point of the BZ (the 
Dirac point of the germanene/silicene), which is now 
above the Fermi level. The gaps that opened for the 
Ger/MoS 2 and Sil/MoS 2 superlattices are 24 and 7 meV, 
respectively (the sizes of the gaps could be well under- 
estimated). Since the electronic states around K point 
are almost fully contributed from the germanene/silicene 
layers, the gaps that opened for the superlattices are due 
to the interactions between the germanene/silicene 
layers only. In other words, the formation of the small- 
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Figure 3 Band structures of free-standing, (a) Germanene calculated with a 4x4 supercell, (b) MoS 2 monolayer calculated with a 5 x 5 
supercell, and (c) silicene calculated with a 4x4 supercell. (d, e) The band structures of Ger/MoS 2 and Sil/MoS 2 superlattices, respectively. The 
contributions from the germanene/silicene and MoS 2 layers to the band structures of the superlattices are shown with blue and green dots, 
respectively. The detailed band structures in the vicinity of the opened band gap are inserted. Red dashed lines represent the Fermi level. 



sized band gaps at the K point is due to the symmetry 
breaking within the germanene/silicene layers caused by 
the introduction of the MoS 2 sheets in the formation of 
superlattices [43-46]. 

To further explore the bonding nature and the charge 
transfer in the Ger/MoS 2 and Sil/MoS 2 superlattices, the 
contour plots of the charge density differences (Api) on 
the planes passing through germanene, silicene, and sul- 
fur layers (in the x-y plane) are shown in Figure 4a,b,c, 
d. The deformation charge density Api is defined as 

^Pi(r) =p(^)-Z!^atom( 7 -^)' where p(j) re P re ' 

sents the total charge density of the superlattice and 

^p atom (V -R^J is the superposition of atomic charge 

ft 

densities. The deformation charge density shown in 
Figure 4a,b,c,d exhibited that the formation of the Ger/ 
MoS 2 and Sil/MoS 2 superlattices did not distort sig- 
nificantly the charge densities of germanene, silicene, 
or sulfur layers, when compared with the deformation 
charge density in the free-standing germanene, silicene 
layers, or sulfur layers in the MoS 2 sheets (not shown). 



Figure 4e,f shows the contour plots of Api on the 
planes perpendicular to the atomic layers and passing 
through Mo-S, Ge-Ge, or Si-Si bonds in the Ger/MoS 2 
and Sil/MoS 2 superlattices. As in the case of isolated 
germanene/silicene or MoS 2 monolayer (not pre- 
sented), the atomic bonding within each atomic layer 
in both the superlattices are mainly covalent bonds. 
Moreover, shown in Figure 4g,h, we also present the 
charge density differences (Ap 2 ) of the same planes as 

in Figure 4e,f. The Ap 2 is defined as Ap 2 = p - 

Pslab( Ger / Sil )-Pslab( MoS 2)> where p(j) > Pslab(Ger/Sil), 

and /? s iab(MoS 2 ) are the charge densities of the super- 
lattice, the germanene/silicene, and the MoS 2 slabs, re- 
spectively. In the calculation of p s i ab (Ger/Sil) and p slah 
(M0S2), we employ the same supercell that is used for 
the superlattice. For calculating the p s iab(Ger/Sil), the 
M0S2 slabs in the superlattice are removed and the 
charge densities of the germanene/silicene slabs are 
then calculated including a structure relaxation. For 
calculating p s i a b(MoS 2 ), the germanene/silicene layers 
are then removed. Such a Ap 2 can clearly demonstrate 
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the charge transfer between the stacking layers in the 
superlattices. Figure 4g,h indicates that the charge 
transfer happened mainly within the germanene/sili- 
cene and the MoS 2 layers (intra-layer transfer), as well 
as in some parts of the intermediate regions between 
the germanene/silicene and MoS 2 layers (inter-layer 
transfer). This is somewhat different from the gra- 
phene/MoS 2 superlattice, where the charge transfer 
from the graphene sheet to the intermediate region 
between the graphene and MoS 2 layers is much more 
significantly visible [6]. Such charge redistributions in the 
Ger/MoS 2 and Sil/MoS 2 systems, shown in Figure 4, indi- 
cate that the interactions between some parts of the stack- 
ing atomic layers are relatively strong, suggesting much 
more than just the van der Waals interactions between the 
stacking sheets. 



Conclusions 

In summary, the first principles calculations based on 
density functional theory including van der Waals cor- 
rections have been carried out to study the structural 
and electronic properties of superlattices composed 
of germanene/silicene and MoS 2 monolayer. Due to 
the relatively weak interactions between the stacking 
layers, the distortions of the geometry of germanene, 
silicene and MoS 2 layers in the superlattices are all 
relatively small. Unlike the free-standing germanene or 
silicene which is a semimetal and the MoS 2 monolayer 
which is a semiconductor, both the Ger/MoS 2 and 
Sil/MoS 2 superlattices exhibit metallic electronic prop- 
erties. Due to symmetry breaking, small band gaps are 
opened up at the K point of the BZ for both the super- 
lattices. Charge transfer happened mainly within the 
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Figure 4 Contour plots of the deformation charge density (Ap n and Ap 2 ). (a, b) kp ] on the planes passing through germanene and sulfur 
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planes perpendicular to the atomic layers and passing through Mo-S, Ge-Ge, or Si-Si bonds in the superlattices. (g, h) Charge density differences 
(Ap 2 ) of the same planes as those in (e) and (f). The green/blue, purple, and yellow balls represent Ge/Si, Mo, and S atoms, respectively. Orange 
and blue lines correspond to Ap > 0 and Ap < 0, respectively. 
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germanene/silicene and the MoS 2 layers (intra-layer 
charge transfer), as well as in some parts of the inter- 
mediate regions between the germanene/silicene and 
MoS 2 layers (inter-layer charge transfer). Such charge 
redistributions indicate that the interactions between 
some parts of the stacking layers are relatively strong, 
suggesting more than just the van der Waals interac- 
tions between the stacking sheets. 
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